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SECTION  I 
INTRODUCTION 


Finite  element  method  is  now  well  established  as  one  of  the  foremost 
numerical  techniques  for  analyzing  linear  elastic  crack  problems. 

Among  other  advantages  associated  with  the  use  of  finite  element  method 
is  the  fact  that  complicated  geometries  and  physical  loading  situations 
can  be  quite  accurately  modeled.  Therefore,  the  method  is  useful  in 
analyzing  standard  as  well  as  nonstandard  crack  growth  and  fracture 
specimens  which  may  be  used  in  the  laboratory. 

A  disadvantage  of  the  method  is  that  with  even  a  slight  change  in  the 
geometry  of  a  given  specimen,  such  as  the  crack  length,  a  new  finite 
element  mesh  is  required.  Since  finite  element  mesh  generation  is 
usually  a  time  consuming  process,  parametric  studies  involving  any  change 
in  geometry  are  tedious.  Using  general  purpose  large  finite  element 
computer  codes,  such  as  NASTRAN  or  MARC,  for  such  problems  is  an  inefficient 
and  expensive  process.  However,  if  special  purpose  finite  element  codes 
developed  primarily  for  parametric  studies  of  test  specimens  are  used, 
the  method  can  be  employed  economically  and  with  improved  efficiency. 

The  present  report  provides  a  user’s  guide  for  a  special  purpose  finite 
element  code  developed  primarily  for  two  dimensional  linear  elastic 
analysis  of  test  specimens.  User's  guides  for  some  supporting  computer 
programs,  such  as  mesh  generators  for  commonly  used  test  specimens  and 
for  plotting  a  mesh,  are  also  Included.  Precise  instructions  and  a  number 
of  illustrative  examples  are  provided  to  clarify  the  data  preparation  for 
executing  the  program. 

It  is  expected  that  a  user  with  only  a  basic  knowledge  of  theory  of  elasticity 
and  some  acquaintance  with  FORTRAN  language  should  be  able  to  effectively 
employ  the  procedures  included  in  this  report. 


SECTION  II 
PROGRAM  HIGHLIGHTS 


1.  ASSUMPTIONS 

1)  Material  is  a  homogeneous,  isotropic,  and  linearly  elastic  solid. 

2)  Conditions  of  plane  stress  or  plane  strain  exist. 

3)  Compared  to  the  dimensions  of  the  solid,  displacements  are  small,  and 
assumptions  of  linear  theory  of  elasticity  are  valid. 

2.  METHOD  OF  SOLUTION 


The  numerical  technique  used  by  the  program  is  based  on  displacement  finite 
element  formulation. 

3.  ELEMENT  GEOMETRY 

The  element  used  in  the  present  program  is  an  eight  node  isoparametric 
quadrilateral  3hown  in  Fig.  1.  Some  special  and  degenerate  forms  of 
this  element  are  shown  in  Fig.  2.  Element  nodes  are  always  numbered 
anticlockwise  starting  at  any  corner  node.  The  four  corner  nodes  are 
number  (I),  (I  +  1),  (1+2),  and  (1+3).  The  first  intermediate  node 
(I  +  4)  is  always  between  corner  nodes  (I)  and  (I  +  1) .  Other  intermediate 
nodes  are  then  numbered  anticlockwise  also. 

4.  MESH  GENERATION 


Mesh  generation  can  be  best  explained  by  considering  a  simple  example: 
consider  the  cantilever  beam  shown  in  Fig,  3.  For  simplicity  let  us  decide 
to  use  only  three  elements.  Fig.  4  shows  some  possible  discretizations . 

The  mesh  of  Fig.  4a  is  most  appropriate  of  all  the  other  meshes  shown.  As 
a  rule  of  thumb  if  the  geometry  of  the  problem  so  allows,  curved,  irregular 
and  degenerate  shapes  should  be  avoided. 
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The  node  numbering  should  be  such  that  the  difference  between  any  two  node 
numbers  belonging  to  the  same  element  is  kept  as  smai:  as  possible.  This 
difference  in  Fig.  Aa  is  7  (*8  -  1).  The  above  practice  helps  in  reducing 
computational  time. 

Elements  may  be  numbered  in  any  arbitrary  way.  However  it  is  a  good  practice 
to  be  systematic  as  in  Fig.  Aa.  Size  of  each  element,  and  the  total  number 
of  elements  in  the  mesh  are  decided  depending  upon  the  nature  of  a  given 
problem.  In  general  a  finer  mesh  would  provide  more  accurate  results  than  a 
coarser  mesh. 

Mesh  generation  can  be  accomplished  either  by  actually  drawing  the  elements 
on  a  graph  paper  or  by  writing  a  mesh  generation  program.  Usage  of  mesh 
generation  programs  for  some  specimen  geometries  is  described  in  Appendix  I. 

5.  TOPOLOGY 


Complete  topology  of  the  mesh  is  described  by  individual  element  connectiv¬ 
ities  and  X  and  Y  coordinates  of  the  nodal  points  with  respect  to  any 
conveniently  chosen  X,  Y  coordinate  system,  as  shown  in  the  following 
example. 


For 

the  mesh 

shown  in  Fig. 

Aa 

the  topology 

can  be 

described 

as  follows. 

a. 

Nodal  Coordinates 

NODE 

X 

Y 

NODE  '  X 

Y 

NODE 

_X 

Y 

1 

A.O 

0.5 

7 

3.50 

1.0 

13 

1.50 

1.5 

2 

A.O 

1.0 

8 

3.50 

1.5 

1A 

1.25 

0.5 

3 

A.O 

1.5 

9 

2.50 

0.5 

15 

1.25 

1.5 

A 

3.75 

0.5 

10 

2.50 

1.5 

16 

1.00 

0.5 

5 

3.75 

1.5 

11 

1.50 

0.5 

17 

1.00 

1.0 

6 

3.50 

0.5 

12 

1.50 

1.0 

18 

1.00 

1.5 
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b.  Connectivities 


Following  the  numbering  convention  described  in  3,  the  element  connectivites 
are  as  follows. 


ELEM 

I 

1+1 

14-2 

1+3 

1+4 

1+5 

1+6 

1+7 

1 

1 

3 

8 

6 

2 

5 

7 

4 

2 

8 

13 

11 

6 

10 

12 

9 

7 

3 

16 

11 

13 

18 

14 

12 

15 

17 

Note 

that  I 

can  be  chosen 

to  be 

any  corner 

node. 

6.  PLOTTING  THE  MESH 

After  generating  a  mesh  and  punching  out  the  nodal  coordinates  and  connec¬ 
tivities  on  data  cards  it  is  often  desirable  to  check  for  any  topographical 
errors.  Perhaps  the  easiest  way  is  to  plot  the  final  mesh.  Usage  of  a 
mesh  plotting  computer  program  developed  for  this  purpose  is  described  in 
Appendix  II. 

7.  SPECIFYING  NODAL  DISPLACEMENT  AND  FORCE  BOUNDARY  CONDITION 

For  simplicity  the  program  (SRL01A) is  designed  to  accept  only  zero  displace¬ 
ment  and  non-zero  force  boundary  conditions. 

For  the  beam  of  Fig.  3  and  using  the  mesh  of  Fig.  4a,  we  can  specify  zero 
X  and  Y  displacement  components  at  nodes  1,  2,  and  3.  Force  of  negative  P 
can  be  specified  in  Y  direction  at  node  18.  The  procedure  for  preparing 
boundary  condition  data  is  described  in  Section  III. 

8.  SPECIFYING  NON-ZERO  STRESS  AND  PRESSURE  ON  ELEMENT  SIDES 


In  problems  where  non-zero  stress  or  pressure  is  to  be  specified  (such  as 
in  crack  problems  with  crack  face  pressure  distribution) ,  the  pressure  or 


stress  distribution  over  any  side  of  the  element  is  given  in  polynomial  form. 
Constant,  linear,  quadratic,  and  cubic  distributions  are  permitted.  The 
method  for  specifying  these  distributions  is  discussed  in  Section  III. 

9.  ENFORCING  CRACK  TIP  STRESS  SINGULARITY 

The  elastic  stress  singularity  at  the  crack  tip  is  included  by  surrounding 

the  crack  tip  by  elements  of  the  form  shown  in  Fig.  2e.  The  collapsed 

nodes  (I  +  1),  (I  +  2),  and  (I  +  5)  are  placed  at  the  crack  tip  and  the 

relative  displacements  among  these  nodes  are  forced  to  be  zero.  This  is 

accomplished  either  by  specifying  the  X  and  Y  displacements  of  the  above 

three  nodes  to  be  zero,  or  simply  by  assigning  a  single  number  to  all 

three  nodes  in  the  nodal  numbering  of  the  mesh.  Details  of  this  method  can 
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be  found  in  the  paper  by  Barsoum.  The  use  of  the  above  procedure  is 
illustrated  in  Section  IV. 


SECTION  III 
DATA  PREPARATION 

1.  USER'S  GUIDE  FOR  SRL01A 

In  this  section  the  method  for  inputting  topology,  material  properties, 
and  boundary  conditions  data  for  the  solution  of  a  given  problem  is 
described.  Each  data  set  refers  to  a  specific  information  required 
by  SRL01A  for  execution. 


a.  Data  Set  1  (sizing  card) ,  Format  (1214) ,  Number  of  cards  »  1 


Columns 

Variable 

Definition 

1-4 

NPRTYP 

Problem  Type  (1  or  2). 

5-8 

NPOINT 

Total  number  of  nodes  in  the  mesh. 

9-12 

NELEM 

Total  number  of  elements  in  the  mesh  . 

13-16 

NBOUN 

Total  number  of  nodes  with  displacement  boundary 
condition. 

17-20 

1 

Total  number  of  nodes  with  concentrated  forces. 

21-24 

1 

Degrees  of  freedom  per  node  (always  give  2) . 

25-28 

NYM 

Number  of  different  materials  in  the  same  problem. 

29-32 

NBAND 

Expected  band  width. 

33-36 

IR 

Number  of  nodes  per  element  (always  give  8) . 

37-40 

NSTRSS 

Stress  computation  indicator  (0  or  1) . 

41-44 

NCPOIN 

Total  number  of  corner  nodes. 

45-48 

NQPTS 

Total  number  of  re-specified  nodes. 

Instructions 


NPRTYP  Give  1  for  Plane  Stress  problem. 

Give  2  for  Plane  Strain  problem. 

NP01NT  Count  and  specify  total  number  of  nodes  in  the  mesh  (corner 
nodes  +  intermediate  nodes). 
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NELEM  Count  and  specify  total  number  of  elements  in  the  mesh. 

NBOUN  Count  and  specify  total  number  of  nodes  which  are  fixed  in  either 

one  or  both  X  and  Y  directions. 

NCONC  Count  and  specify  total  number  of  nodes  which  have  concentrated 
applied  forces  in  either  one  or  both  X  and  Y  directions. 

NYM  In  general  each  element  in  the  mesh  can  be  allowed  to  possess 

different  elastic  properties.  If  so,  specify  the  total  number 
of  materials  to  be  used.  For  uniform  elastic  properties  for  all 
elements  give  NYM  equal  to  1. 

NBAND  This  is  a  sizing  parameter  found  by  the  following  formula. 

NBAND  *  (Max.  difference  between  any  two  nodes  of  an  element 

+  1)  x  2. 

NSTRSS  For  crack  problems  where  stress  intensity  factors  are  of  interest 
give  NSTRSS  =  o.  For  other  problems  where  element  stress  and 
strain  components  are  required  give  NSTRSS  =  1. 

NCPOIN  Total  number  of  corner  nodes  (NPOINT  -  number  of  intermediate  nodes) 
in  the  mesh. 

NQPTS  If  NQPTS  is  given  zero,  the  intermediate  nodes  are  automatically 

placed  at  the  middle  of  straight  line  distance  between  the  adjoining 
corner  nodes.  For  problems  in  which  this  arrangement  needs  to  be 
altered,  such  as  in  crack  problems  using  the  element  of  Fig.  2e, 
the  number  of  those  intermediate  nodes  which  have  to  be  moved 
from  the  mid-position  should  be  counted  and  specified  as  NQPTS. 

b.  Data  Set  2  (Elastic  Property  Cards),  Format  (2F20.5),  Number  of 


Columns 


Variable 


Definition 


Young's  Modulus 
Poisson's  Ratio 


See  a.  Section  III. 


c.  Data  Set  3  (Boundary  Conditions) .  Format  (414) .  Number  of  Cards  ■  NBOUN 


Columns 

Variable 

Definition 

1-4 

NF(1) 

Node  with  displacement  boundary  condition 

5-8 

NB(I,1) 

Zero  X-displacement  Indicator 

9-12 

NB(I,2) 

Zero  Y-displacement  indicator 

Instructions 

NF(I)  Give  the  number  of  a  node  whose  X  or  Y  or  both  displacements  are 
fixed. 

NB(I,1)  Give  1  if  X-displacement  of  node  NF(I)  is  fixed.  Give  0  if 
X-displacement  of  node  NF(I)  is  not  fixed. 

NB(I,2)  Same  as  NB(I,1)  but  for  Y-displacement  of  node  NF(I). 


See  a.  Section  III. 


d.  Data  Set  4  (Concentrated  Loads).  Format  (I4.2F20-1Q) .  Number  of 
Cards  =  NCONC* 


Columns 

Variable 

Definition 

1-4 

5-24 

25-44 

NP 

U(NPx2-l) 

U(NPx2) 

Node  number  of  point  with  cone,  load 
Value  of  X-component  of  applied  load 
Value  of  Y-component  of  applied  load 

See  a.  Section  III. 


Data  Set  S  (Corner  Point  Coordinates).  Format  (I4,2F10*5),  Number 


of  Cards  -  NCPOIN 


Columns 


Variable 

Definitions 

I 

** 

Corner  node  number 

X(I) 

X-coordinate  of  I 

Y(I) 

Y-coordinate  of  I 

See  a.  Section  III. 

!* 

For  example  in  Fig.  4a  points  1,  3,  6,  8,  11,  13,  16,  and  18  are  corner 
nodes (NCPOIN  -  8). 


Data  Set  6  (Element  Connectivities),  Format  (9I4,F10*5,I4) ,  Number 


of  Cards  ”  NELEM 


Columns 


1-4 

5-8 

9-12 

13-16 

17-20 

21-24 

25-28 

29-32 

33-36 

37-46 

47-50 


Variable  I  Definitions 


First  corner  node 

Second  corner  node 

Third  corner  node 

Fourth  corner  node 

1+4  First  Intermediate  node 

1+5  Second  intermediate  node 

1+6  Third  intermediate  node 

1+7  Fourth  intermediate  node 

NEP  Material  property  number 

THICK  Element  Thickness 

NTYPEL  Applied  pressure/stress  indicator 

(0  or  1) 


Instructions 


NEP  If  a  single  material  is  used  (NYM  *  1  in  Data  Set  1),  NEP  is 

1  for  all  elements. 

NTYPEL  If  any  side  of  an  element  (see  Fig.  1)  is  subjected  to  applied 

pressure,  for  that  element  give  NTYPEL  as  1.  For  other  elements 
NTYPEL  is  given  0. 

*See  a.  Section  III. 

Follow  the  numbering  convention  of  Part  3.  (Section  II) 


g.  Data  Set  7  (Crack  Tip  Node) ,  Format  (14) ,  Number  of  Cards  ■  1 


Columns 

Variable 

Definition 

1-4 

NTIP 

Node  number  of  a  node  located  at 
the  crack  tip 

If  there  is  no  crack  give  zero. 


h.  Data  Set  8  (Coordinate  Modification),  Format  (I4.2F10-5).  Number  of 
Cards  «  NQPTS* 


Columns 

Variable 

Definition 

1-4 

I 

Node  number  of  the  point  whose 

coordinates  are  to  be  modified 

5-14 

X(I) 

X  -  coordinate  of  I 

15-24 

Y(I) 

Y  -  coordinate  of  I 

* 

See  a. 


Section  III. 
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L. 


Data  Set  9  (Pressure  Polynomial),  Format  (14),  Number  of  Cards  ~  M 

Format  (8F10»5),  Number  of  Cards  *» 


Card  1 


M 


Columns 

Variable 

Definition 

1-4 

NS1DE 

Side  number  of  element  on  which  pressure 
is  applied  (See  Fig.  1) 

Card  2 


Columns 


Variable 


Definition 


** 


A1 

A2 

21-30 

a3 

31-40 

a4 

41-50 

»1 

51-60 

b2 

61-70 

b3 

71-80 

®4 

First  coefficient  of  traction  polynomial 

Second  coefficient  of  traction  polynomial 
Third  coefficient  of  traction  polynomial 

Fourth  coefficient  of  Ty  traction  polynomial 
First  coefficient  of  Tx  traction  polynomial 
Second  coefficient  of  Tx  traction  polynomial 
Third  coefficient  of  Tx  traction  polynomial 
Fourth  coefficient  of  Tx  traction  polynomial 


This  is  an  optional  data  set.  Use  only  if  one  or  more  elements  have 
NTYPEL  =  1  (See  f..  Section  III).  If  there  are  no  such  elements  omit  this 
data  set  altogether. 

* 

M  is  the  number  of  elements  with  NTYPEL  »  1.  For  each  such  element  there 
are  two  cards.  The  set  of  Card  1  and  Card  2  has  to  be  given  for  each 
such  element.  Give  M  such  sets.  The  order  of  the  sets  is  the  same  as  the 
order  of  element  cards  in  Data  Set  6  which  have  NTYPEL  =  1. 

Ty  =  Ai  +  A2X  +  A3X2  +  A4X3 
Tx  -  Bi  +  B2Y  +  B3Y2  +  B4Y3 

(Continued  at  top  of  next  page) 
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T  and  T  are  the  x-  and  y-  traction  components  respectively,  obtained 
x  y 

by  taking  the  vector  dot  product  between  the  stress  components  and  the 
outward  unit  normal  on  the  surface. 


j.  Data  Set  10  (Kt  Computation  Indicator) .  Format  (14) ,  Number  of  Cards 


If  NTIP-0  (see  g.  Section  III),  this  card  is  omitted, 
k.  Data  Set  11  (Kl  Computation  Point),  Format  (14),  Number  of  Cards  -  1* 


Variable 

Definition 

NKIC 

Nodal  point  on  the  crack  edge  whose 
displacement  is  to  be  used  for  Ki 
calculation.  NKIC  should  be  approximately 
A/ 50.0  distance  away  from  NTIP,  where  A 
denotes  the  crack  length 

If  NTIP-0) (see  g.  Section  III),  this  card  is  omitted. 

2.  LIMITATIONS 

For  using  the  stored  version  of  SRL01A  the  following  limitations  on  the 
values  of  various  variables  shduld  be  observed. 

0  <  NPOINT  <  340* 

0  <  NELEM  <  100 
0  <  NBOUN  <  50 
0  <  NBAND  <  100 
0  <  NYM  <  5 

* 

For  definitions  of  these  variables  see  a.  Section  III. 
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SECTION  IV 
SOLUTION  PROCEDURE 


1.  PROGRAM  EXECUTION 

The  necessary  control  cards  for  executing  the  program  with  a  single  data 
set  (one  problem)  or  multiple  data  sets  (many  problems)  are  as  follows. 

First  Card  (Job  Card),  start  at  column  1. 

AAA, CM300000 , T100 , 10100 , STANY .  M760328 

AAA  m  Job  identification  number  (eg.  user's  initials). 

T100  •  Execution  time  limit  in  seconds.  This  may  be  changed  (eg.  T60) 

10100  -  Input-Output  time  limit  in  seconds.  This  may  be  changed. 

M7 60328  -Example  of  user's  code.  Use  appropriate  account  number. 

Note:  Put  a  period,  then  a  blank  after  STANY. 

Second  Card  (attach  card),  start  at  column  1. 

ATTACH, LGO.SRLOIA,  ID-M760328. 

SRL01A  ■  Program  name.  For  other  programs  described  in  appendices  the 

program  name  is  changed  but  other  parameters  remain  identical.  For 
executing  SRL01A  finite  element  program  use  the  above  card  exactly 
the  way  it  is  shown. 

Note  the  period  after  the  ID  value. 

Third  Card  (load-go  card),  start  at  column  1. 

LGO , 
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Fourth  Card  (end-of-record) ,  this  card  is  available  at  the  computer  center. 


CARD  SEQUENCE 

The  necessary  control  cards  and  data  cards  have  to  be  arranged  in  the 
following  sequence. 

Job  card 
Attach  Card 


LGO.  N 


LGO.  J 


As  many  as  number  of  problems  to  be  solved. 


End  of  record  card 
Data  for  first  problem 
End  of  record  card 
Data  for  second  problem 


] 

} 


End  of  record  card 
Data  for  the  last  problem 
End  of  job  card. 


} 
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SECTION  V 

NUMERICAL  EXAMPLES 

The  purpose  of  providing  solved  numerical  examples  using  the  programs 
discussed  in  the  present  manual  is  to  further  clarify  the  data  preparation 
and  program  execution  described  in  previous  sections. 

1.  EXAMPLE  1:  SINGLE  NOTCH  SEMI-CIRCULAR  SPECIMEN 

In  order  to  explain  all  the  procedural  details ,  let  us  start  with  the 
problem  of  a  semi-circular  ring  containing  a  radial  crack  shown  in  Fig.  6. 
Due  to  axial  symmetry  about  the  X-axis  only  one  half  of  the  geometry  is 
needed  for  analysis  (Fig.  6b).  To  account  for  the  other  half,  we  need  to 
specify  Y-displacement  (V)  to  be  zero  on  the  boundary  indicated  in  Fig.  6b. 

The  purpose  is  to  obtain  nodal  displacements  and  stress  intensity  factors 
(Ki)  for  various  values  of  crack  length  (A) ,  inner  radius  (Ri) ,  and  outer 
radius  (R0) .  We  proceed  in  the  following  steps: 

a.  Step  1:  Mesh  Generation 

As  explained  in  Section  2,  mesh  generation  can  be  accomplished  either  by 
hand  or  by  using  the  computer.  Let  us  use  SRLC  program  of  Appendix  I 
to  generate  the  mesh.  The  two  data  cards  used  were: 

Card  1 :  Blank 

Card  2:  0-5  1-0  2-0 

The  numbers  in  card  2  represent  A,  R^,  R0,  respectively. 
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The  output  of  SRLC  in  the  form  of  punched  deck,  and  in  print,  was  as 
follows : 


•5U003  l.uUQOQ  2.00600 
1US  52  17  1  2  1  70  8  0  of  5 

lOCdoCOG.uCuOC  • 3U030 

10?  0  1 

10  6  0  1 

107  0  1 

108  b  1 

109  0  1 

114  0  1 

123  6  1 

12d  U  1 

137  0  1 

142  0  1 

151  0  1 

156  w  1 

165  0  1 

170  0  1 

179  0  1 

184  0  1 

185  1  1 

13  1.  OOPCatiOOOO 


1 

•  bjllvO 

l.oOOCO 

3 

•  Jot-1 

1.1125b 

5 

«  uuGUU 

1.22500 

7 

• UQOOQ 

1.50000 

9 

. tocoo 

1. 7750  C 

11 

•  0  o  0  0  0 

1.88750 

13 

•  ii  U  U  bl> 

2.  000  OL 

21 

-.382fa8 

•92386 

23 

-.42574 

1.02782 

25 

-•40679 

1.13175 

27 

-•  5741,3 

1. 38582 

29 

-•  67926 

1.63989 

31 

-.72231 

1.74382 

33 

-.76537 

1.84776 

41 

-.70711 

.70711 

43 

-. 7  o 6o6 

•  7  8Lub 

45 

— . 66621 

•86621 

47 

— 1. 061 66 

1.06066 

49 

-1.26511 

1.29511 

51 

-1.334oo 

1. 33b66 

53 

-1.41421 

1.41421 

61 

-.92368 

.  362  66 

63 

-1. C  2762 

.42574 

65 

-1.13175 

.45879 
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67 

•1.38582 

.57403 

69 

-1.53969 

•67926 

71 

-1.74382 

.72231 

73 

-1.84776 

•76537 

81 

-• 990 wO 

•  11C  00 

83 

-1.11250 

.11000 

85 

-1.22500 

.11000 

87 

-1.50000 

•  11800 

89 

-1.77500 

.11000 

91 

-1. 80760 

.11000 

93 

-1* 98uu  J 

• llbCO 

108 

-1*  QOGwO 

0. OOU00 

102 

-1.11250 

a. Goooa 

106 

-1.225ofl 

G. 00009 

105 

-1.775uO 

0.  OCiooO 

107 

-1.88750 

0.00000 

109 

—  2. 0  j o  j 0 

0.00000 

115 

-1. 337?o 

0. UOuuO 

117 

-1. 337»0 

.  0  9  0  o  3 

119 

-1. 5ulo0 

.  o90l)Q 

121 

- 1.  30  2  90 

. 09u00 

123 

-1. 33290 

o  •  0 00 o  9 

129 

—1.35375 

u.UuoUO 

131 

-1.393/6 

• o75oQ 

133 

-1.  90  0  0  0 

.  O?5C0 

135 

“i. 30625 

• 0  75oG 

137 

-1.  oJ6<;5 

0 .  oOoQO 

143 

-1.43000 

a.ooouo 

145 

-1. 490  O0 

.05000 

1h7 

-1. 9W&00 

.  o5 j  Ou 

149 

-1. 33000 

.  05000 

151 

-1. 990  OU 

0 . 0  0  0  J  0 

157 

-1. 48o oo 

0. OuouQ 

169 

-1*  980 wo 

. 02000 

161 

-1. 9U0  J  0 

•  02000 

163 

- 1.  52  0  oU 

.  020  UO 

165 

-1.520ou 

0. OGOuQ 

1/1 

-1.99uoO 

0. 00000 

173 

-1.  99293 

. 0U7U7 

17» 

-1.90000 

.01000 

177 

-1. 30  7  j  7 

. 00707 

179 

-1. 510  oO 

o .  o  o  o  o  o 

165 

-1. 5u0o0 

o .  J  0  U  u  0 

3 

23  21 

1 

15  22 

5 

25  23 

3 

16  24 

7 

27  25 

5 

17  26 

9 

C'i  2  7 

f 

18  26 

11 

31  25 

9 

1 3  30 

13 

33  31  11 

20  32 

14 

2 

i 

1. J9300 

15 

4 

1 

1 .  J  0  3  o  u 

lo 

o 

1 

1. JC J03 

17 

8 

1 

1.  J  0  J  J  u 

16 

10 

1 

1.  U  0  0  0  y 

19 

12 

1 

'  1.30  3 OG 

23 

h3 

*1 

21 

35 

42 

34 

22 

1 

1. 

jjOud 

2$ 

49 

40 

23 

3o 

H  4 

36 

24 

X 

1* 

J  J  tf  b  6 

2/ 

*7 

■»9 

29 

37 

4o 

3b 

26 

1 

1. 

J  0  0  6  0 

29 

49 

♦  7 

27 

3o 

48 

37 

26 

1 

1. 

0  j  306 

31 

91 

>♦9 

29 

39 

50 

36 

3u 

1 

1. 

0  J  3  u  u 

33 

93 

51 

31 

46 

62 

39 

32 

1 

1. 

J  j  J  3  u 

ti 

u3 

OX 

4X 

65 

62 

94 

42 

1 

X. 

j  j  3  *  0 

*♦5 

03 

o3 

43 

56 

64 

65 

44 

1 

1. 

bO  3.0  D 

47 

ot 

09 

45 

57 

66 

5b 

46 

1 

1. 

JbObC 

49 

o9 

o7 

47 

96 

66 

57 

46 

1 

1. 

JO  Ob  3 

9l 

7  x 

3  9 

**9 

3  9 

70 

96 

50 

1 

1. 

J  j  3  j  0 

53 

73 

71 

51 

6  6 

7  2 

59 

52 

X 

1. 

J  *  3  0  6 

63 

83 

61 

61 

76 

82 

74 

62 

1 

1. 

JO  0  0b 

65 

09 

30 

03 

76 

64 

75 

64 

1 

1. 

JuObO 

o7 

67 

69 

69 

77 

66 

76 

06 

1 

1. 

J  J  16  b 

69 

69 

o7 

b7 

78 

66 

7  7 

68 

1 

1. 

jjota 

71 

9i 

69 

69 

79 

96 

78 

7  u 

1 

1. 

J  J  J  J  b 

73 

93 

91 

7l 

80 

92 

79 

72 

1 

1. 

03  60b 

81 

63 

162 

let 

62 

9o 

lbl 

94 

1 

1. 

6  u  3  J  0 

63 

65 

164 

162 

8* 

96 

103 

95 

1 

1. 

jOIuO 

69 

91 

1*7 

1  65 

96 

96 

lub 

97 

A. 

1. 

0  j  0  3  4* 

91 

93 

169 

167 

32 

99 

Uo 

98 

1 

1. 

J  J  3  b  0 

65 

117 

119 

164 

111 

1x6 

x  10 

96 

1 

X. 

30063 

69 

6/ 

119 

1 17 

6b 

112 

116 

111 

1 

1. 

36  6  3  0 

67 

69 

1X1 

119 

08 

113 

120 

112 

1 

1. 

J  J  3  0  0 

121 

o9 

1 6  9 

123 

113 

9/ 

114 

122 

1 

1. 

03  300 

117 

131 

129 

119 

125 

130 

124 

116 

1 

1. 

3  b  0  bk 

117 

119 

133 

131 

118 

126 

132 

125 

1 

1. 

j  j  J  0  3 

119 

121 

l3o 

1  33 

12U 

127 

134 

126 

1 

1. 

OOOuC 

135 

121 

123 

137 

127 

122 

128 

13b 

1 

1. 

OuOOO 

131 

149 

143 

129 

139 

144 

138 

130 

X 

X. 

JOOjO 

131 

133 

147 

145 

132 

146 

146 

139 

1 

1. 

06  ObO 

133 

139 

145 

147 

134 

141 

146 

140 

1 

X. 

JuJuC 

149 

135 

137 

151 

141 

136 

1  42 

150 

1 

1. 

j  3  0  0  0 

145 

199 

19? 

143 

193 

156 

152 

144 

1 

1. 

Ow  ObO 

147 

149 

lo3 

lbl 

146 

155 

162 

154 

1 

X. 

00303 

145 

147 

lol 

159 

146 

154 

160 

153 

X 

1. 

0  3  3  06 

163 

149 

lol 

165 

195 

150 

156 

164 

1 

1. 

Ori  3  0  C 

159 

173 

171 

157 

167 

172 

166 

156 

1 

1. 

00300 

159 

lol 

175 

173 

166 

168 

174 

167 

1 

1. 

06000 

161 

163 

177 

176 

162 

169 

176 

168 

1 

1. 

33  300 

177 

163 

1 99 

179 

169 

164 

176 

178 

1 

1. 

03360 

173 

169 

169 

171 

161 

165 

180 

172 

1 

1. 

00006 

175 

165 

185 

173 

182 

165 

161 

174 

1 

1. 

00063 

177 

165 

165 

175 

183 

169 

162 

176 

1 

1. 

30063 

179 

165 

160 

161 

182 

163 

164 
1 

171 

169  185  177 

-1. 49756 
-1*49823 
-1.96600 
-1.56177 
-1.9J250 

184  189 

6* 66600 
•00177 
.00256 
.60177 
6.06006 

163 

178 

1 

1. 

06  000 
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The  first  line  in  the  output  gives  the  crack  length  and  the  inner  and 
outer  radii  of  the  specimen.  The  rest  of  the  computer  print-out  corresponds 
to  the  data  sets  described  in  Section  III  as  follows. 


Line  2 

Data 

set 

1 

Line  3 

Data 

set 

2 

Lines  4  to  20 

Data 

set 

3 

Line  21  ' 

Data 

set 

4 

Lines  22  to  88 

Data 

set 

5 

Lines  89  to  140 

Data 

set 

6 

Line  141 

Data 

set 

7 

Lines  142  to  146 

Data 

set 

8 

Line  147 

Data 

set 

10 

Line  148 

Data 

set 

11 

Note  that  data  set  9  is  missing  because  no  element  side  is  subjected  to 
applied  pressure. 


Mesh  Plotting 


From  the  punched  data  deck  of  step  1,  data  sets  5  and  6  (coordinates  and 
connectivities)  were  temporarily  removed  for  use  in  the  plotting  program 
SRL11.  In  accordance  with  the  description  provided  in  Appendix  II, 
the  following  data  was  supplied  to  SRLll. 


6  185  52 

[Nodal  coordinates  of  data  set  5  above] 

[Blank  card] 

[Connectivities  of  data  set  6  above] 

1  20  1  C-SPEC 

0.0  2.0  -2.0  -1.0  3.0  3.0  1.0  1.0 

[Blank  card] 

End -of -record. 
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1 


9 

1 


By  executing  SRL11,  the  plot  shown  In  Fig.  7  was  obtained.  In  Fig.  8 
details  of  the  fine  mesh  near  the  crack-tip  are  shown.  Fig.  8  was  obtained 
by  supplying  SRL11  with  only  those  element  connectivities  which  were  to 
be  plotted  and  by  increasing  the  scales  (XSCALE,  YSCALE) .  The  number  of 
elements  (52)  in  the  first  data  card  corresponds  to  the  number  of  elements 
which  are  to  be  plotted. 

c.  Step  3:  Executing  SRL01A 

Data  sets  5  and  6  were  placed  back  in  their  proper  sequence  in  the  punched 
data  deck  of  Step  1  after  removing  the  blank  card  which  was  inserted 
behind  data  set  5  for  plotting  purpose.  SRL01A  was  then  executed  by 
using  the  following  cards. 

Job  card 

ATTACH ,LGO , SRL01A , ID-M7  60  3 28 . 

LGO. 

End -of -record 
[Data  deck  of  Step  1] 

End-of-job. 

The  output  of  SRL01A  was  as  follows. 
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The  Kj  value  for  the  problem  is  computed  by  the  program  and  is  given  at 
the  end  of  the  output.  Displacement  values  at  various  locations,  such 
as  crack  mouth  opening  displacement  (node  100)  and  load  point  displacement 
(node  13)  can  be  read  directly  by  identifying  the  proper  node  number  on 
the  mesh.  It  should  be  noted  that  the  actual  crack  mouth  opening  displace¬ 
ment  and  total  load  line  displacement  in  the  present  case  will  be  twice 
the  Y-displacement  of  nodes  100  and  13,  respectively.  The  results  are: 

Kj  *  21-58609  psi/In  ,  Load  Point  Displacement  =  233-74620  x  10"?  in. 
Crack  Mouth  Opening  Displacement  =  75.1878  x  10~7  £n 

e.  Remarks 

AAA, CM700G , T100 , 10100 , STANY .  M - 

ATTACH,  TAPES,  SRLC1,  ID-M760328 
ATTACH, LGO,  SRLC.ID-M760328 . 

LGO. 

End-of- record 

3  (Total  number  of  cases  minus  1) 

0.2  1.0  2.0 

0.3  1.0  2.0 

0.4  1.0  4.0 

0.5  1.0  4.0 

End-of -job . 

The  above  set-up  will  produce  data  for  crack,  lengths  0.2,  0.3,  0.4,  and 
0.5  for  Rj  *  1.0  and  R<>  ■  2.0.  The  job  set-up  for  multiple  runs  using 
SRL01A  is  described  in  Section  IV. 

Note  that  the  last  four  elements  are  as  shown  in  Fig.  2e.  Nodes  180,  181, 
182,  183,  and  184  are  located  at  quarter  points  of  the  element  sides.  This 
necessitates  the  specification  of  the  coordinates  of  these  points  in  data 
set  8  (Section  III).  The  purpose  of  this  particular  procedure  is  to  model 
the  proper  stress  singularity  or  the  crack  tip  as  mentioned  in  Section  II. 
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2.  EXAMPLE  2:  DOUBLE  NOTCH  RING  TENSION  SPECIMEN 

Consider  a  double  notch  circular  ring  (Fig.  9)  under  diametral  tension. 

Due  to  symmetry  about  X  and  Y  axes  only  one  quarter  of  the  ring  needs  to 
be  considered  for  analysis.  In  fact  by  changing  only  the  boundary  conditions 
the  same  mesh  as  of  Example  1  (Fig.  7)  can  be  used.  This  Is  accomplished 
simply  by  changing  the  data  sets  1,  3,  and  4  of  Example  1.  The  following 
computer  print-out  shows  the  changes  in  the  above  data  sets. 
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Note  that  in  data  set  1  the  number  of  boundary  conditions  (NBOUN)  has 
been  changed  from  17  to  30.  In  data  set  3,  thirteen  more  boundary 
conditions  have  been  added  and  the  boundary  condition  of  the  crack  tip  has 
been  changed.  In  data  set  4,  the  load  Is  now  applied  at  node  1.  Also, 
since  the  applied  load  is  1*0  lbs,  the  resulting  values  for  displacements 
and  Ki  should  correspond  to  a  2*0  lbs  load  on  the  actual  ring. 

The  rest  of  the  solution  procedure  remains  identical  to  that  of  Example  1. 
The  computer  print-out  of  the  results  obtained  by  SRL01A  is  shown  on  the 
following  pages. 
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As  in  Example  1,  the  Kj  value  corresponding  to  the  applied  load  (2*0  lbs) 
is  given  directly.  Nodes  1  and  100  correspond  to  the  load  point  and  the 
crack  mouth,  respectively.  Accordingly,  we  have  for  1  lb.  load: 

Kx  -  3 • 0068  psi/In 

load  point  displacement  *  22*34021  x  10-7  in. 

crack  mouth  opening  displacement  «  9*511695  x  10~7  in. 

3.  EXAMPLE  3:  DOUBLE  NOTCH  RING  COMPRESSION  SPECIMEN  (CRACK  LINE  PRESSURE) 

The  geometry  of  a  Double  Notch  Compression  (DNC)  Specimen  Is  shown  in 
Fig.  10a.  Let  us  choose  A  -  0*5  in,  Rf  -  1*0  in. and  R<>  “  2*0  in. so  that 
the  mesh  of  examples  1  and  2  may  again  be  utilized.  For  this  case  however, 
instead  of  obtaining  a  solution  for  a  concentrated  load  directly,  as  in 
Example  2,  let  us  use  the  crack  line  pressure  concept;  the  purpose  being 
to  illustrate  the  procedure  for  applying  distributed  loads. 

In  Fig.  10b  cjyy  is  the  crack  line  pressure  obtained  by  analyzing  an  unflawed 

ring  (no  cracks)  under  the  given  loads  P.  In  general,  crack  line  stress 

Oyy  may  be  obtained  either  by  solving  the  analytical  elasticity  problem  or 

by  performing  finite  element  stress  analysis  of  the  unflawed  ring.  The 

analytical  solution,  if  available,  is  of  course  preferable.  For  the  case 

under  consideration,  the  stress  distribution  on  the  crack  line  was  found  by 

3 

taking  the  negative  of  the  analytical  stress  distribution  obtained  for  the 
unflawed  ring  when  P  *  1.  The  node  numbers,  nodal  coordinates  and  ir.CTyy  of 
the  points  lying  on  the  crack  line  are  listed  below. 

Node _ 

100  (crack  mouth) 

101 
102 

103 

104 


Node 

X-coordinate (in) 

ir.Oyy(ps£) 

110 

-1.28125 

-2.68151 

115 

-1.3375 

-1.94734 

124 

-1.36563 

-1.61820 

129 

-1.39375 

-1.30913 

138 

-1.42188 

-1.01697 

143 

-1.45 

-0.739051 

152 

-1.465 

-0.595854 

157 

-1.48 

-0.455750 

166 

-1.485 

-0.409688 

171 

-1.49 

-0.363927 

180 

-1.4975 

-0.295833 

185  (crack  tip) 

-1.5 

-0.273276 

It  is  noticed  that  nodes  100,  101,  and  102  belong  to  the  same  element  whose 
connectivity  is:  81  83  102  100  82  95  101  94  (see  connectivity  data  of 
Example  2).  In  accordance  with  the  convention  shown  in  Fig.  1,  the  nodes 
100,  101,  and  102  form  side  3  of  the  element.  In  a  similar  fashion  we 
can  identify  the  group  of  nodes  (102,  103,  104),  (104,  110,  115), 

(115,  124,  129),  etc.  being  side  3  of  respective  elements  to  which  they 
belong . 


Next,  for  each  element  side,  such  as  (100,  101,  102),  we  have 
polynomial  traction  distribution  of  the  form: 

Ty  *  A1  +  A2X  +  A3X2  +  A$X^  =  Uyy  ny  ■  Oyy 

0 


Tx  "  B1  +  ®2  +  B3Y2  +  B4Y3 


where  n^  is  the  only  nonzero  component  of  the  outward  unit  normal.  For  the 
present  case  Bj  -  B2  *  B3  ■  B4  ■  0,  Assuming  a  quadratic  distribution 
of  Ty,  A4  *  0.  Then,  for  each  group  of  nodes,  such  as  (100,  101,  102), 
the  coefficients  A^,  A3  are  found  since  the  coordinates  and  pressure 
distribution  at  each  node  is  known.  For  the  present  case  this  was 
accomplished  by  writing  a  small  computer  program  to  solve  three  linear 
algebraic  simultaneous  equations.  The  result  was  as  follows. 


NODE  GROUP 


2 


HDI 

100,101,102 

51-03588 

80-16584 

32-35651 

102,103,104 

28-82489 

39.99275 

14.19187 

104,110,115 

17.87724 

22.03494 

6.82783 

115,124,129 

12.93540 

14.57838 

4.01532 

129,138,143 

10.66238 

11.31630 

2.84494 

143,152,157 

9.28895 

9.41632 

2.18785 

157,166,171 

8.69650 

8.61401 

1.91623 

171,180,185 

8.42768 

8.25336 

1.79527 

Using  the  displacement  boundary  conditions  shown  in  Fig.  10b,  the  SRL01A 
data  for  the  problem  can  now  be  prepared.  Since  the  mesh  and  displacement 
boundary  conditions  remain  identical  to  those  of  Example  2,  only  the  changes 
are  described  in  the  following. 


a.  Data  set  1:  NCONC  *  0  (no  concentrated  loads).  . 

b.  Data  set  2:  Same  as  Example  2. 

c.  Data  set  3:  Same  as  Example  2. 

d.  Data  set  4:  Not  present. 

e.  Data  set  5:  Same  as  Example  2. 

f.  Data  set  6:  NTYPEL  =  1  for  those  elements  along  the  crack  surface. 

g.  Data  set  7:  Same  as  Example  2. 

h.  Data  set  8:  Same  as  Example  2. 

i.  Data  set  9: 
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51.03588 

a 

80.16584 

32.35651 

0.0 

0.0 

0.0 

0.0 

0.0 

j 

28.82489 

'X 

39.99275 

14.19187 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

17.87724 

1 

22.03494 

6.82783 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

12.93540 

14.57838 

4.01532 

0.0 

0.0 

0.0 

0.0 

0.0 

3 


10.66238 

O 

11.3163 

2.84494 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

9.28895 

9.41632 

2.18785 

0.0 

0.0 

0.0 

0.0 

0.0 

3 

8.6965 

8.61401 

1.91623 

0.0 

0.0 

0.0 

0.0 

0.0 

J 

8.42768 

8.25336 

1.79527 

0.0 

0.0 

0.0 

0.0 

0.0 

Data  set 

10:  Same  as 

Example  2, 

Data  set 

11 :  Same  as 

Example  2, 

• 

The  data  image  produced  by  SRL01A,  and  the  finite  element  solution  were 
as  follows. 
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THE  CRACK  TIP  IS  AT  NOOl  165 

”*  STHLSj  intensity  factors 

USING  THt  Y-OISPL  OF  NOPt  171  Kl=  2. *0167 


The  above  results  show  that  Kx  for  the  problem  is  2*60167  psiVTn.  The 
same  problem  of  double  notch  ring  compression  specimen  when  solved  by 
applying  concentrated  load  P  directly  (Fig.  10a)  gave  stress  intensity 
factor  value  of  2*59355  psi/in.  The  0.3%  difference  in  the  two  values 
could  be  attributed  to  round-off  error  in  the  calculation  of  nodal  stresses 
from  elasticity  solution  and  in  the  evaluation  of  pressure  polynomial 
constants.  Further,  the  assumption  of  quadratic  pressure  distribution 
over  each  element  side  may  have  contributed  to  the  difference.  An  interesting 
exercise  would  be  to  try  constant,  linear,  and  cubic  pressure  polynomials. 

4.  EXAMPLE  4:  COMPACT  TYPE  SPECIMENS 

One  of  the  advantages  of  using  programmed  mesh  generation  along  with  the 
finite  element  program  is  that  changes  in  geometry  pose  no  difficulty. 

For  a  compact  tension  test  specimen  (Fig.  11)  of  standard  geometry 
(H/W  -  0.6)  and  for  WOL  specimen  (H/W  -  0.486)  SRL01A  results  can  be 
found  in  reference  4.  For  the  present  example  we  choose  the  following 
dimensions: 

H  -  1.2  in. 

W  =  2.4  in. 

WX  -  3.05  in. 

S  *  0.0938  in. 

F  »  1.417  in. 

E*  -  0.55  in. 

R  -  0.25  in. 

9  »  40.0  deg. 

A  *  1.2  in. 
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The  mesh  generation  program  SRLCMP  (Appendix  I)  was  executed  using  the 
following  cards: 

Job  Card 

ATTACH , LGO , SRLCMP , ID-M760328 . 

ATTACH, TAPE8 , SRLCM, CY-2 , ID=M760328 . 

LGO. 

End-of-record  card 

1*2  0.55  1*417  1*2  0*0938  40*0  2*4  3*05  0*25  0*55 

End-of-job  Card 

The  last  entry  on  the  data  cards  is  the  distance  n  -  n  (Fig.  11).  This 
distance  should  always  be  kept  less  than  2  *  E’  and  larger  than  S.  The 
location  of  point  n  can  be  chosen  to  correspond  with  the  actual  COD 
measurement  location  in  an  experimental  set-up. 

The  use  of  plotting  program  SRL11  provided  the  mesh  plots  shown  in 
Fig.  12.  Due  to  symmetry  about  the  crack  axis,  only  the  upper  half  of 
the  geometry  is  used  in  the  analysis.  The  displacement  and  force 
boundary  conditions  are  shown  in  Fig.  13.  The  crack  tip  is  fixed  in  both 
X  and  Y  directions  to  eliminate  rigid  body  motion. 

The  data  image  produced  by  SRL01A  was  as  follows: 
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The  results  of  the  analysis  are  summarized  below: 

Kx  =  6-59341  psi/In 
^m-m  =  66-43846  x  10-7  In 
6 1-1  -  43-76998  x  10*7  in 
6p_p  -  48-98702  x  10"7  in 
6g_g  -  46-30266  x  10"7  in 

The  points  m,  1,  p,  and  g  are  indicated  in  Fig.  11  and  these  correspond 
to  nodes  318,  261,  277,  and  281  on  the  mesh  respectively.  The  displace¬ 
ment  values  shown  above  were  obtained  by  multiplying  the  V-displacements 
of  the  corresponding  points  by  two. 

5.  EXAMPLE  5:  THREE  POINT  BEND  SPECIMENS 


The  geometry  of  a  three  point  bend  specimen  is  shown  in  Fig.  14.  The 
data  for  this  problem  was  generated  by  using  the  program  of  subsection  4 
of  Appendix  I,  and  the  following  dimensions: 

L  -  4-25  in-  ,  Pf  =  0-0625  in. 

S  -  4.0  in.  ,  THETA  =  40.0  deg. 


S  *  4.0  in. 

W  =  1.0  in . 

G'  =  0.375  in. 
N  =0.25  in. 


G  =  0.2  in. 

A  =  0.5  in. 

H  *  0.125  in. 


B  *  1.0  in 


The  data  image  produced  by  SRL01A  was  as  follows: 
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362  366  316  316  361  367  315  368"  1  1  1*66666 
-M6_299  .312  316  299  366  313  367  1  1  1.66666 
312  298  296  316  366  297  365  311  1  1  1461666 
-316  -336 -332  316  322  331  323  317  1  1  1.66666 
316  316  328  336  315  321  329  322  1  1  1.69666 

JL14  3A2L  126  328  313  326  327  321  -  1  _  1 . 1.66666 

326  312  310  326  326  311  319  325  1  1  1*66666 
-336  338  366  332  336  339  337  331  1  1  1*86666 
328  338  360  336  335  339  336  329  1  1  1966606 
J26  336  366  328  336  339  335  327  1  1  1*66666 
326  336  366  326  333  339  336  325  1  1  1968686 
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NOCE  199 
NOCE  209 
NOCE  291 
NOCE  202 
NODE  203 
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NOCE  212 
NOCE  213 
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NOCE  217 
NOCE  216 
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In  this  case  loads  of  0.5  lbs.  are  applied  at  points  1  and  p  (Fig.  14), 
and  the  crack  tip  is  fixed  in  both  directions.  Therefore  the  total  load 
point  displacement  is  the  sum  of  the  magnitudes  of  vertical  displacements 
of  points  p  and  1.  In  Fig.  14,  G  represents  the  location  of  clip  gage 
to  measure  the  crack  opening  displacement.  The  result  of  the  analysis 
gave  the  following  values  of  stress  intensity  factor,  load  point  displace¬ 
ment  and  the  crack  opening  displacement: 

Kl  =  10.6892  psi/in 

Load  point  displacement  =  58.01097  x  10“  ^  in. 

Crack  opening  displacement  =  16.02759  x  10”?  in. 

On  the  mesh  points  p  and  1  correspond  to  nodes  30  and  290,  respectively. 
Crack  opening  displacement  is  twice  the  Y  displacement  of  node  229. 

6.  EXAMPLE  6:  SINGLE  NOTCH  RING  SPECIMEN 

Single  notch  specimens  in  the  shape  of  circular  rings  are  shown  in  Fig.  15. 
The  program  of  subsection  2  of  Appendix  I  generates  the  data  for  single 
notch  ring  compression  specimen  (Fig.  15a).  The  data  for  the  tension 
case  is  generated  by  changing  the  loading  condition. 

Choosing  A  =*  0.5  in.,  =  1.0  in  and  RQ  =  2.0  in.,  the  mesh  generation 
program  was  used  to  produce  the  required  data  for  SRL01A.  The  data  image 
produced  by  SRL01A,  which  corresponds  to  the  mesh  of  Fig.  16  was  as 
follows . 
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The  execution  of  SRL01A  resulted  In  the  following  values  for  stress 
intensity  factor,  crack  opening  displacement,  and  load  point  displacement: 

Kx  =  1.8940  psi/In 

COD  =  2  x  (Y-displacement  of  node  180)  *  6.45932  x  10”^  in. 

6l,p  =  sum  of  the  absolute  value  of  X  -  displacements  of  node  189 
and  node  13  =  15.31517  x  iq-7  in. 
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SECTION  VI 
CONCLUSION 


The  report  contains  user's  manuals  for  the  following  programs: 

1.  SRL01A:  A  two  dimensional  linear  elastic  finite  element 

analysis  code  for  crack  problems. 

2.  SRL11:  A  mesh  plotting  program. 

3.  SRLC:  A  mesh  generation  program  for  C-shaped  and  double¬ 

notch  ring  shaped  specimens. 

4.  SRLRNG:  A  mesh  generation  program  for  single-notch  ring 

shaped  specimens. 

5.  SRLCMP:  A  mesh  generation  program  for  compact  tension  type 

specimen  geometries. 

6.  SRLBND:  A  mesh  generation  program  for  bend  specimen  geometry. 

Although  the  mesh  generation  programs  have  been  carefully  tested,  the 
user  should  not  expect  these  programs  to  produce  acceptable  data  for 
arbitrarily  chosen  geometrical  dimensions  for  any  given  specimen  shape. 

The  programs  are  however  expected  to  generate  proper  data  for  standard 
geometries  with  minor  modifications. 

It  is  highly  recommended  that  a  new  user  try  to  solve  one  or  more  of  the 
illustrative  examples  presented  in  Section  IV  before  attempting  to  solve 
a  new  problem. 
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APPENDIX  I 

MESH  GENERATING  PROGRAMS 


1.  C-SHAPED  SPECIMEN 

The  finite  element  data  for  SRL01A  to  analyze  a  C-shaped  specimen  may  be 
generated  by  executing  the  following  program: 

Job  Card 

ATTACH, LGO, SRLC, ID=M760328 . 

ATTACH, TAPE8 , SRLC1 , ID=M760328 . 

LGO. 

End-of-record 

N  (Total  number  of  data  sets  to  be  generated  minus  1) 

A  RI  RO  (Crack  length,  inner  radius,  outer  radius) 


A  RI  RO 
End-of-job. 

Note:  a)  Format  for  N  is  (14) 

b)  Format  for  A,  RI,  RO  is  (3F10.5) 

c)  The  mesh  may  also  be  used  to  analyze  double  notch  circular 
ring  specimens  (see  examples  2  and  3). 

2.  SINGLE-NOTCH  RING  SPECIMEN 


Data  for  a  single  notch  ring  specimen  (Fig.  15a)  can  be  generated  by 
using  the  following  program. 
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Job  Card 

ATT  aCH  ,  LGO ,  SRLRNG ,  ID-117  60328 . 

ATTACH , TAPE8 , SRLRG , CY=2 , ID-M76032  8 . 

LGO. 

End-of-record 

A  RI  RO  (crack  length,  inner  radius,  outer  radius) 

End-of-job. 

Note:  a)  Format  for  the  data  card  is  (3F10.5) 

b)  For  multiple  runs  put  as  many  LGO.  cards  as  number  of  cases, 
and  place  an  end-of-record  card  in  front  of  each  data  card. 

c)  The  mesh  may  also  be  used  for  a  single  notch  ring  tension 
specimen  (Fig.  15b)  simply  by  changing  the  location  of  the 
applied  load  from  nodes  13  and  189  to  node  81. 

3.  COMPACT  TYPE  SPECIMEN 


The  data  to  analyze  a  compact  specimen  of  given  geometry  may  be  generated 
by  using  the  following  program. 

Job  Card 

ATTACH , LGO , SRLCMP , ID-M760328 . 

ATTACH, TAPE8 , SRLCM, CY-2 . ID-M760328. 

LGO. 

End-of-record 

A  E  F  H  S  THETA  W  W1  R  GS 
End-of-job 

Note:  a)  Format  for  the  data  card  (10F7.4) 

b)  GS  Is  the  distance  n-n  in  Fig.  11. 

c)  Procedure  for  multiple  runs  is  the  same  as  for  single  notch 
ring  specimen  of  the  previous  section. 
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4.  THREE  POINT  BEND  SPECIMEN 


Note 


Job  Card 

ATTACH , LGO , SRLBND , ID=M7 60328. 

ATTACH, TAPE8,SRLB,CY=2,ID=M760328. 

LGO. 

End-of-record 

L  S  W  G'  N  P'  THETA  G  A  H  B 
End-of-j ob 

a)  Format  for  the  data  card  (11F7.4)  for  definition  of  symbols 
see  Fig.  14. 

b)  For  multiple  runs  follow  the  procedure  described  in  Section 
3  for  single  notch  ring  specimen. 
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APPENDIX  II 
PLOTTING  PROGRAM 

1.  USER'S  GUIDE  FOR  MESH  PLOTTING  PROGRAM  (SRL11) 

a.  Data  Set  1,  Format  (314),  Number  of  cards  =  1 


— ~  — 

Columns 

Variable/Constant 

Definition 

4 

6 

Output  unit 

5-8 

NPOIN 

4e 

Total  number  of  nodes 

9-12 

NELEM 

* 

Total  number  of  elements 

f  * 

f  See  a.  Section  III. 

i 

i 


b.  Data  Set  2  ,  Form  (I4,2F10*5),  Number  of  cards  =  NCPOIN 


Columns 

Variable 

Definition 

1-4 

I 

Node  number  of  corner  point 

5-14 

X(I) 

Y-coordinate  of  I 

15-24 

Y(I) 

Y-coordinate  of  I 

Use  the  same  data  cards  as  5.  Section  II. 
** 

See  a.  Section  III. 


c.  Data  Set  3,  Format  (not  applicable).  Number  of  cards  =  1 
(Blank  Data  Card) 
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1.  Data  Set  4,  Format  (814),  Number  of  cards  =  NELEM 


Columns 

Variable 

Definition 

1-4 

I 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

5-8 

I  +  1 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

9-12 

1  +  2 

Element  connectivities.  Use  Data 
Set  6  of  f.,  Section  III. 

13-16 

1  +  3 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

17-20 

1  +  4 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

21-24 

1  +  5 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III, 

25-28 

1  +  6 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

29-32 

1  +  7 

Element  connectivities.  Use  Data 
Set  6  of  f..  Section  III. 

See  a.  Section  III. 


Data  Set  5,  Format  (315, A10) ,  Number  of  cards  =  1 


T 


T 


Columns 


Variable/Constant 


Definition 


AD-AO07  MO  SYSTEMS  RESEARCH  LABS  INC  DAYTON  OHIO  F/6  9/2 

TWO  DIMENSIONAL  LINEAR  ELASTIC  ANALYSIS  OF  FRACTURE  SPECIMENS  U — ETC(U) 
FEB  00  J  AHMAD  F33615-79-C-502S 

UNCLASSIFIED  AFWAL-TR-80-4000  NL 


2  >  2 

■ 

END 

9-80 

OTIC 

I 


Instructions 


L  Give  1  if  elements  and  nodes  are  not  to  be  numbered. 

Give  2  if  only  elements  are  to  be  numbered. 

Give  3  if  nodes  and  elements  are  to  be  numbered. 

Give  4  if  only  nodes  are  to  be  numbered. 

LABEL  Give  ten  character  title  (including  spaces,  such  as 

RING  SPMN1). 

NOTE:  Only  the  corner  nodes  of  the  mesh  are  numbered  by  the  program. 


f .  Data  Set  6,  Format  (8F10‘5),  Number  of  cards  =  1 


Columns 

Variable/Constant 

Definition 

i-10 

YMIN 

Minimum  of  Y-coordinates  of  all  nodes 

11-20 

YMAX 

Maximum  of  Y-coordinates  of  all  nodes 

21-30 

XMIN 

Minimum  of  X-coordinate  of  all  nodes 

31-40 

XMAX 

Maximum  of  X-coordinates  of  all  nodes 

41-50 

YSCALE 

* 

Scale  on  the  Y  axis 

51-60 

XSCALE 

* 

Scale  on  the  X  axis 

61-70 

1*0 

71-80 

1*0 

It  is  recommended  to  use  SCALE  ■  YSCALE  ■  Plot  dimension  desired  divided 
by  actual  mesh  dimension.  Plot  dimension  desired  is  limited  by  the  size 
of  the  paper  available  at  CALCOMP  plotter.  Usually,  Y  dimension  of  the 
plot  »  (YMAX  -  YMIN)  *  YSCALE  should  not  exceed  8  inches. 


g.  Data  Set  7.  Format  (not  applicable).  Number  of  cards  ■  1 

(Blank  Data  Card) 
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2.  SAMPLE  JOB  SET-UP  FOR  THE  BEAM  PROBLEM  OF  SECTION  II.  5. 

AAA , CM150000 , T100 , IOIOO , STANY .  M760328 
ATTACH , CCAUX , CCAUX, ID-LIBRARY , SN-ASD 
LIBRARY, CCAUX. 

ATTACH, LGO , SRL11 , ID-M7  60328 . 

LGO. 

End-of-record  card. 


6 

18 

3 

1 

4.0 

0.5 

3 

4.0 

1.5 

6 

3.5 

0.5 

8 

3.5 

1.5 

11 

1.5 

0.5 

13 

1.5 

1.5 

16 

0.5 

0.5 

18 

0.5 

1.5 

.ANK 

CARD 

1 

3 

8 

8 

13 

11 

16 

11 

13 

1  12 

0.5 

1.5 

1.0 

BLANK  CARD 
End-of-data  card 


BEAM 

.0  2.0 
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Figure  1.  Eight  Noded  Quatrilateral  Element. 


Figure  4.  Some  Possible  Discretizations 


Figure  12  (a).  Mesh  for  Compact  Specimen. 

Upper  Half  of  Spec  Inert. 


i 


Figure  12  (b) . 


Mesh  for  Compact  Specimen. 
Details  of  Crack-Tip  Region. 
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YOUNGS  MODULUS  -  I07PSI 
POISSON'S  RATIO  »  0.3 
THICKNESS  »  1.0  IN. 


Figure  13.  Upper  Half  of  Compact  Specimen.  (Schematic) 
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(a) 


Figure  IS.  Single  Notch  Ring  Specimens 

(a)  Tension 

(b)  Compression 
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